International Europhysics Conference on HEP 



PROCEEDINGS 



Steps towards the QCD calculation for e+e 3 jets 
at NNLO 



Sven Moch, Peter Uwer 

Institut fiir Theoretische Teilchenphysik, Universitdt Karlsruhe, 
76128 Karlsruhe, Germany 
, E-mail: noch@particle.uni-karlsruhe.de, iiwerQparticle .uni-karlsruhe .de| 



O 

CN ' Stefan Weinzierl* 



Dipartimento di Fisica, Universita di Parma, INFN Gruppo Gollegato di Parma, 



Q . 43100 Parma, Italy 



E-mail: 



^ m 

stefanw@fis.unipr.it 

cn: iXl 

Abstract: High precision analyses of experimental data for e+e annihillation, such as 

^ \ determination of jet rates or event shape observables, call for complete next-to-next-to- | ■ 

I leading order (NNLO) perturbative QCD predictions. In this talk, we discuss the various ' 

■"sj" ' ingredients entering the calculation of the NNLO corrections in e^e~ annihillation. CD 

o: ^ 

I ■ 1. Introduction o 

(— I \ Perturbation theory is a powerful tool for precise theoretical predictions on the outcome of 

high-energy experiments. The state-of-the-art is the transition to fully differential next-to- 

■ next-to-leading order (NNLO) calculations for jet physics. The most prominent processes 

■ where a complete NNLO calculation is desirable are Bhabha scattering, — > 2 jets and 



e"'"e~ — > 3 jets. A NNLO calculation of e'^e" — > 3 partons is expected to reduce the 
theoretical uncertainty in the extraction of down to 1% 0. Furthermore, it allows to 
model the jet structure more accurately and should improve the knowledge on the interplay 
between perturbative and power corrections. 

In the past years there has been tremondous progress towards this goal: The relevant 
master integrals for — > 2 jets and e'^e~ 3 jets have been calculated, either with the 
help of a Mellin-Barnes representation Q or by solving differential equations |^] . Numerical 
results for these master integrals serve as a useful cross-check |Q, ^ . With these results the 
two-loop amplitudes for Bhabha scattering Q, for ^ 2 jets [^, ^, for pp — > 77 [^] and 



for light-by-light scattering |10| have been calculated. 
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However, the two-loop amplitudes are only one part of the story. Loop amplitudes in 
a massless theory like QCD lead to infrared singularities. These divergences cancel against 
corresponding singularities resulting from amplitudes with additional unresolved partons. 
Presently, the singular behaviour of double- unresolved tree amplitudes [11, ^] and of one- 



loop amplitudes with one unresolved parton is understood |13, 111], as well as the structure 



of the poles of two-loop amplitudes up to 1/e [15|. But despite this progress there is still 
a fair amount of work to be done: 

• The two-loop amplitude for e'^e~ — > 3 jets remains to be calculated, although recently 



there has been reported progress fl^. Due to multiple scales, the amplitude will 
depend on the ratio of two kinematical invariants, for example S12/S123 and 823/5123- 
In contrast, the two-loop amplitudes calculated so far depend only on one ratio of 
invariants, say s/t. 

• A method to cancel IR divergences has to be set up. This requires the extension of 
the subtraction or slicing method to NNLO. In particular analytic integrations over 
unresolved regions in phase space have to be carried out. 

• The final numerical computer program requires stable and efficient numerical meth- 
ods, in particular for the Monte Carlo integration of the double unresolved contribu- 
tions. 

In this talk we focus on the first two problems. 
2. Two-loop amplitudes 

Let us briefiy review how the calculation of the two-loop amplitudes for the single scale 
problems has been performed. Starting from the Feynman diagrams, which yield two-loop 
tensor integrals, i.e. integrals with powers of the loop momentum in the numerator, one 
performs the tensor reduction. Using a Schwinger parametrization these tensor integrals 
can be related to scalar integrals with higher powers of the propagators and different values 



of the dimension D [17]. In a second step, integration-by-parts ]18], Lorentz-invariance 



identities ]y] or form- factor relations [19\ can be used to eliminate propagators and to 
reduce the integrals to a few basic topologies. At this stage, the basic topologies occur in 
various dimensions and with various powers of the propagators. The last step relates 
these basic topologies to the (much smaller) set of master integrals, which are in most cases 
just the basic topologies for D = 4 — 2e and all = 1. Again, integral relations derived 
with integration-by-parts or form-factor relations are used for this step. Looking at the 
computational cost, one finds that the first two steps are rather easy to perform, whereas 
the third step requires to solve large systems of equations. In particular in the presence of 
multiple scales this becomes rather involved. 

We therefore would like to advocate a different approach, where basic topologies with 
arbitrary powers of the propagators and arbitrary dimensions are evaluated directly ]pO|] . 
As an advantage, this approach keeps the size of the intermediate expressions much smaller 
and allows for a very efficient implementation of the reduction scheme. The starting point 
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is the observation that these integrals can be written as nested sums involving Gamma 
functions. For example, a basic integral occuring in e^e~ 3 jets is the so-called C- 



topology, which can be written as [20| 



r(2m - 2£ - ;vi235)r(l + 1^1235 - 2m + 2e)T{2m - 2e - z^2345)r(l + z^2345 - 2m + 2g) 

r(i/i)r(i^2)r(z^3)r(z^4)r(z^5)r(3m -3e - 2^12345) 



X 



— X 



r(m - £ - U5)T{m - e - 1^23) ,_ .2m-2e-vu3A5 

r(2m - 2£ - Z.235) ^ .^o^fe)^i'^2! 

r(ii + z^3)r(i2 + V2)^{h + ^2 - 2m + 2£ + z/i2345)r(ii + i2-m + e + 1/235) 
r(ii + 1 - 2m + 2£ + z/i235)r(i2 + 1 - 2m + 2£ + V2M5)T{ii + ^2 + ^^23) 

,2m-2e-i/i235 



r(ii + 2m - 2£ - z/i25)r(«2 + z^2)r(^i + ^2 + z^4)r(ii + i2 + - e - z^i) 
^ r(ii + 1 + 2m - 2£ - z/i235)r(i2 + 1 - 2m + 2£ + z/2345)r(ii + ^2 + 2m - 2£ - i/15) 

_ 2m-2£-1^2345 

r(ii + i/3)r(i2 + 2m - 2£ - z/345)r(n + «2 + z^i)r(n + «2 + "t- - e - 2^4) 

^ r(ii + 1 - 2m + 2£ + z/i235)r(i2 + 1 + 2m - 2£ - z/2345)r(ii + i2 + 2m - 2£ - 2^45) 

2m-2£-i.i235 ^2m-2£-i.2345 r(n + 2m - 2£ - Z/l25)r(i2 + 2m - 2£ - t/345) 



r(ii + 1 + 2m - 2£ - z/i235)r(i2 + 1 + 2m - 2£ - z/2345 
r(ii + ^2 + 2m - 2£ - i/235)r(ii + ^2 + 3m - 3£ - 2^12345 



r(ii + i2 + 4m - 4£ - 1/12345 - 1^5) 



2.1) 



where we set xi = (— si2)/(— S123) and X2 = (— S23)/(— si23)- Here, the result for the 
integral holds for arbitrary dimensions D = 2m — 2£ and any (not necessarily integer) 
power Vi of the propagators. 

The Gamma functions can be expanded systematically in £, thus allowing to solve 
these nested sums algorithmically to any given order in e and to express the result in the 
basis of Z-sums, defined by 

Z{n;mi,...,mk;xi,...,Xk) = V -4ir----4r- (2-2) 

— ^1 

n>n>j2>...>«fc>0 

Important special cases of this definition for Z-sums are multiple polylogarithms (in which 



we express our final results) ||2l| - [p3 | 



Umk,...,mi{xk, .■.,xi) = Z(oo;mi,...,mfc;xi,...,Xfc) (2.3) 
and Euler-Zagier sums (which occur in the expansion of Gamma functions) 

^mi,...,mfe("') = Z(n;mi,...,mfc;l,...,l). (2.4) 
Z-sums form an algebra. In fact, Z-sums can be considered as generalizations of Euler- 



Zagier sums or harmonic sums. Many algorithms known for the latter |24] carry over 
to Z-sums. The usefulness of the Z-sums lies in the fact, that they interpolate between 
multiple polylogarithms and Euler-Zagier sums, the interpolation being compatible with 
the algebra structure. 
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3. Subtraction method 

The NNLO cross section receives contributions from double unresolved tree amplitudes, 
one- loop amplitudes with one unresolved parton and two-loop amplitudes: 

7 (1) _ / Atree* ^l-loop Al-loop* Atree\ 

da^^^ = (^A'r^ * Al-^""P + Al-^""P* A'^^^ + A\^^°°^* ^„#„, (3.1) 



'n.+l"'Pn+l) 



where Qn is the jet-defining function for n partons. Taken separately, each part gives a 
divergent contribution. Only the sum of all contributions is infrared finite. This problem 
already occurs in NLO calculations and can be solved by adding and subtracting a suitable 



chosen term. For example the NLO cross section for (n -|- 1) partons is written as |25] 

I datl^ + I da'il, = I [da'^l, - da^,) + J + da^^,) . (3.2) 

Here da^_^_l acts as a local counterterm for single unresolved configurations to da^^2 ™ ^ 
dimensions and is integrable over a one-parton subspace. 

A similar subtraction scheme is needed for the NNLO cross section. The more com- 
plicated part is the subtraction term to da^^2 for double unresolved configurations. In 
addition, the one-loop corrections require a subtraction term dan'^°°^ which approx- 

imates in D dimensions the one-loop corrections when one parton becomes unresolved. We 
shortly outline how to obtain the subtraction terms for dan'^°°^ . We first decompose the 
amplitudes into primitive amplitudes (e.g. with a fixed cyclic ordering and a definitive 
routing of the fermion lines). Primitive one-loop amplitudes factorize in the collinear limit 
where one parton becomes unresolved as ||l3| , [ij] 

Atree* A-loop . .1-loop* Atree 

Sing*^^'^ ( ^^^^^ * + Ai'^°°P * At[^^^ + Singi-'°°P ^^^^^ * A^;[^^ (3.3) 

The singular function Sing*'''^^ is already known from NLO calculations and poses no prob- 
lem. The function Sing^~'°°^ gives a new contribution. However, it is relatively simple to 
write down an appropriate subtraction term for each primitive structure. As an example 
we consider the splitting g — > qq. For the primitive part where the fermion entering the 
loop turns right we have as subtraction term 



_gpU ^S^"" 

2piPj 



fi'^ Y ( 1 3 1 1 



Here S^"^ denotes the spin correlation tensor. This subtraction term has to integrated over 
a one-parton phase space. For the more complicated integrals we can rely on the technique 
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of nested sums |26[. The integrated counter-part reads: 



-P2 



11 /23 1 \ 1 /I 2 23 118\ 1 



'3e3 Vl8 + 3 7e2 ^ 18^ ^ 27 J e 

7. 1 3 23 2 118 1 2 4075\ ] ^, , 

We note that in general for one-loop amplitudes with one unresolved parton the soft radi- 
ation pattern is more complicated as compared to a NLO calculation. Generally, a simple 
dipole picture with emitter and spectator will not be sufficient. 



4. Outlook and conclusions 

The NNLO calculation for e^e~ — > 3 jets is a challenging project. We have reported on 
recent progress in the evaluation of two-loop amplitudes with multiple scales using the novel 
approach of nested sums. All algorithms for the solution of the nested sums are suitable 
for implementation in computer algebra systems like GiNaC ||2^ or FORM p^ . We have 
outlined the road to the cancellation of infrared divergences. This involves the analytic 
integration over the single and double unresolved partonic phase space and we believe the 
technique of nested sums to be of great value here as well. We have not addressed any 
issues concerning the final numerical integration of the fully differential NNLO cross 
section by means of Monte Carlo techniques. 
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